Dynamics of coherences in the interacting double-dot Aharonov-Bohm interferometer: Exact 

numerical simulations 
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We study the real time dynamics of electron coherence in a double quantum dot two-terminal Aharonov- 
Bohm geometry, taking into account repulsion effects between the dots' electrons. The system is simulated 
by extending a numerically exact path integral method, suitable for treating transport and dissipation in biased 
impurity models [Phys. Rev. B 82, 205323 (2010)]. Numerical simulations at finite interaction strength are 
supported by master equation calculations in two other limits: assuming non-interacting electrons, and working 
in the Coulomb blockade regime. Focusing on the intrinsic coherence dynamics between the double-dot states, 
we find that its temporal characteristics are preserved under weak-to-intermediate inter-dot Coulomb interaction. 
In contrast, in the Coulomb blockade limit, a master equation calculation predicts coherence dynamics and a 
steady-state value which notably deviate from the finite interaction case. 

PACS numbers: 73.63.Kv,03.65.Yz, 73.63.-b, 73.23.Hk 
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I. INTRODUCTION 

Is electron transfer through quantum dot structures phase 
coherent, or incoherent? How do electron-electron and 
electron phonon interactions affect phase-coherent transport? 
From the other direction, what is the role of the interference 
phenomena on many-body effects, such as the formation of 
the Kondo resonance? These questions were addressed in 
numerous experimental and theoretical works, detecting the 
presence of quantum coherence in mesoscale and nanoscale 
objects, using Aharonov-Bohm (AB) interferometry, see for 
example Refs. lfH- [l5ll . In particular, oscillations in the con- 
ductance resonances of an AB interferometer, with either one 
or two quantum dots embedded in its arms, were demonstrated 
in Refs. JHHl], indicating on the presence of quantum coher- 
ence. Interestingly, AB oscillations were also manifested in 
the co-tunneling regime, implying that phase coherence is in- 
volved within such processes 0|. 

Considering the role of electron-electron (e-e) interactions 
in the AB interferometry, a systematic analysis carried out 
in Ref. Jit] has argued that spin flipping channels of the 
transferred electron, the result of e-e repulsion effects, induce 
dephasing. The consequence of this decohering effect was 
the suppression of AB oscillations and the appearance of an 
asymmetry in the resonance peaks. One should note however 
that this study has assumed infinitely strong e-e interactions 
(Coulomb blockade regime) and treated the system perturba- 
tively in the dot-metals coupling stren gth. In other studies, e-e 
repulsion effects were totally ignored IU6I1 . incorporated using 
a mean-field scheme, see for example H, or treat ed p erturba- 
tively using the Green function formalism lfl7i Il8ll . These 
studies, and other theoretical and numerical works JH QIJ, 
have typically considered only the steady-state limit, analyz- 
ing the conductance, a linear response quantity, or the current 
behavior, often in the infinite large bias case JaHl. 

The coherence of electron transfer processes through an AB 
interferometer has been typically identified and characterized 
via conductance oscillations in magnetic fields. However, in 
a double-dot AB structure, a device including two dots, both 



connected to biased metal leads, it is imperative that the rela- 
tive phase between the two dot states (or charge states) should 
similarly convey information on electron coherence and de- 
coherence, as this phase is tangled with the AB phase. In a 
recent work, Tu et al. 12011 have analyzed this intrinsic coher- 
ence dynamics, revealing the effect of phase localization for 
different magnetic fluxes, by studying the real time dynamics 
of the two-dots reduced density matrix. This analysis, based 
on an exact (nonmarkovian) master equation method |[2"lll . left 
out e-e interaction effects all-together. 
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FIG. 1: Scheme of a double-dot AB interferometer. The two dots 
are each represented by a single electronic level. Electron repulsion 
energy is represented by the double arrow. The total magnetic flux is 
denoted by <J>. 



Detailed study of the dynamical role of finite electron- 
electron interactions on the intrinsic coherence behavior in a 
biased double-dot AB interferometer, is the focus of our work. 
The system includes a parallel quantum dot setup for the AB 
interferometer, where (spinless) electrons experience an inter- 
dot repulsion effect. For a schematic representation, see Fig. 
Q] A unified description of the conductance behavior of this 
model, a steady-state property, was given in Ref. 0. Here, 
we focus on the dynamics of the coherences, off diagonal ele- 
ments of the double-dot reduced density matrix. Furthermore, 
we simulate the charge current in the system, assuming dif- 
ferent values for the magnetic flux, at finite bias. Other ef- 
fects considered are the role of finite temperature on the co- 
herence pattern, and the behavior away from the electron-hole 
symmetric point, a regime not considered before in a non- 
perturbative calculation within the AB setup J22I1 . 
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We follow the nonequilibrium real time dynamics of elec- 
tron coherence in this subsystem-bath model (double quantum 
dot-metals) by performing exact numerical simulations, em- 
ploying the recently developed influence functional path inte- 
gral (INFPI) technique 11231 12411 . This method relies on the ob- 
servation that in out-of-equilibrium (and/or finite temperature) 
cases bath correlations have a finite range, allowing for their 
truncation beyond a memory time dictated by the voltage- 
bias and the temperature. Taking advantage of this fact, an 
iterative-deterministic time-evolution scheme has been devel- 
oped where convergence with respect to the memory length 
can in principle be reached. As convergence is facilitated at 
large bias, the method is well suited for the description of 
the real-time dynamics of mesoscale and nanoscale devices 
driven to a steady-state via interaction with biased leads. The 
INFPI approach is complementary to other numerically exact 
methods such as numerical renormalization group techniques 
HHim, real time quantum Monte Carlo simulations l27ll and 
path integral methods l28ll . It offers flexibility in defining the 
impurity object and the metal band structure. The results well 
converge at large voltage bias and/or high temperatures, as we 
show below. 

The principles of the INFPI approach have been detailed 
in Refs. ll23ll24ll . where it has been adopted for investigating 
dissipation effects in the nonequilibrium spin-fermion model, 
and the population and the current dynamics in correlated 
quantum dots, by investigating the single impurity Anderson 
model j2^| and the two-level spinless Anderson dot l30ll . In 
this paper, we further extend this approach, examining the ef- 
fect of a magnetic flux on the intrinsic coherence dynamics. 
Our simulations show that general dynamical characteristics 
of the double-dot coherence are maintained upon the appli- 
cation of inter-dot Coulombic interactions. In particular, the 
characteristic timescale for reaching the steady-state limit, the 
dependence of the coherence on the AB phase factors, and the 
form of the temporal current, similarly evolve for systems at 
zero or finite inter-dot interaction, for finite bias (beyond lin- 
ear response), away from the electron-hole symmetric point, 
at low or high temperatures. We compare our data to (analytic) 
results based on a master equation treatment. This method can 
readily handle the zero e-e interaction case and the counter 
case, the infinite interaction limit. Interestingly, in the latter 
Coulomb blockade case the coherence is expected to evolve 
and sustain values distinctively different from its behavior at 
finite interactions. 

The paper is organized as follows. In Sec. II we describe 
the model and draw the principles of the INFPI technique. In 
Sec. Ill we present numerical results for the coherence dy- 
namics and the charge current, analyzing the role of electron- 
electron interaction. Sec. IV includes analytic results based 
on master equations. Conclusions follow in Sec. V. 



II. MODEL AND NUMERICAL METHOD 

We focus on the symmetric AB setup, with a quantum dot 
(impurity) located at each arm of the interferometer. The dots 
are each connected to two metal leads (referred to as "baths" 



or "reservoirs"), maintained in a biased state. For simplic- 
ity, we neglect the spin degree of freedom and describe each 
quantum dot by a single spinless electronic level. Overall, the 
dots T and '2' are represented by the electronic levels ei, 
and £2, respectively, described by the creation operators d] n , 
(to = 1,2). These levels are coupled in an AB geometry to 
two metal leads (a = L, R) with chemical potentials [i a . For 
a schematic representation see Fig. Q] The total Hamiltonian, 
H, includes the following terms 

H = Eini + e 2 n 2 + Umn 2 + ^ ffcC^c^ 

a,k 

+ [ V LM, m e^— di lCL , k + V R ^ m e^- R c\ k d m 

fc,m=l,2 

+/i.C.] (1) 

Here, k denotes the creation (annihilation) of an electron 
with momentum k in the a lead. We assume identical leads, 
characterized by the same band structure. For the subsystem, 
n m = d\ n d m represents the number operator for the impurity 
level to, U is the charging energy penalty for a simultaneous 
occupancy at the two dots. The AB phase factors, 4> m ,a, are 
acquired by electron waves under a magnetic field perpendic- 
ular to the device plane, 

(f>l,L - 02, L + 4>l,R - 4>2,R = <fi = 27T$/$ . (2) 

Here $ is the magnetic flux enclosed by the ring and <f>o = 
h/e is the flux quantum. In what follows, we adopt the fol- 
lowing gauge, (f) 1>L - (f> 2 ,L = <l>i,R - <t>2,R. = 0/2. Besides 
the phase factors, the coupling strength V a ,k,m are taken as 
real numbers. The hybridization elements are given by 

r aim ,« = 7T V a ,k, m V aXn e i ^"'-^6{e - e fe ). (3) 

k 

We assume that the couplings are identical for the two levels, 
Vk.fe.m = V at k,n, and define the diagonal decay to the a bath 

r Q = 7r]T(V Q , fc , m ) 2 <5(e-e fe ). (4) 

fe 

The total diagonal decay is denoted by V = Tl + Tr. In 
practice, we take T a to be identical at the two ends. Further, 
in what follows we only consider the degenerate situation with 

e = £m- 

In the absence of magnetic fields this model is referred to 
as the "spinless two-level Anderson model". It has been ex- 
tensively studied in the context of molecular electronics, for 
exploring various effects in molecular conduction I3~i1l . and 
in mesoscopic physics, revealing nontrivial effects such as 
pop ulation inversion IU9l [30ll and transmission phase lapses 

MSI. 

Using the INFPI approach, the following observables could 
be followed: the dots' occupation, (n m ) = Tr[dJ„c? m p], the 
coherence, <j\^, = Tr[c?{c?2p], and the total current passing 
through the interferometer. The trace is performed over all 
degrees of freedom, metals and impurity. The charge current 
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presented will be the symmetrized current, (I e ) = Tr[I e p], 
accessed by defining the operator of interest as 



L = VL,k, m c^ L k d m e~ 

k , rn 



(5) 



with 3 denoting the imaginary part. Within INFPI, these ob- 
servables are simulated in the Heisenberg representation as 
we explain below, assuming an initial density matrix p(0) de- 
scribing a nonequilibrium-biased situation. 

We outline now the principles of the INFPI method, allow- 
ing for the exact simulation of transport and dissipation in im- 
purity models I23I l24tl . We begin by reorganizing the Hamil- 
tonian, Eq. (fl]i, as H = Hq + Hi, identifying the nontrivial 
many-body interaction term as 



Hi = U 



nin 2 - - {ni + n 2 ) 



(6) 



H contains the remaining two-body terms, redefining the dot 
energies as Ed, m = e m + U/2. This partitioning allows us 
to utilize the Hubbard-Stratonovich (HS) transformation Isill . 
see Eq. dU below. Formally, the dynamics of a quadratic oper- 
ator, A, either given in terms of the baths (metals) or impurity 
degrees of freedom, can be written as 



(A(t)) = Tr[ P (0)A(t)} 



lim -?-Tr[p(0)e im e xA e- im ].(7) 



Here A a real number, taken to vanish at the end of the calcu- 
lation, p is the total density matrix, and the trace is performed 
over both subsystem and reservoirs degrees of freedom. For 
simplicity, we assume that at the initial time t = the dots 
and the baths are decoupled, p(0) = er(0) <E> pl <8> Pr- The 
baths are prepared in a nonequilibrium biased state p a \ the 
subsystem is described by the (reduced) density matrix a(0). 

We proceed and factorize the time evolution operator, 
_ ^HSt^N ^ f ur t ner utilizing the Trotter decomposi- 



ifft 



tion e 



HSt 



^ e iH 5t/% e iH^t e iH st/%y The many _body term 
H\ can be eliminated by introducing auxiliary Ising variables 
s = ± via the HS transformation IK34I1 . 



±iHiSt 



,H±(s) — -SK±(n 2 -n 1 ) 



(8) 



Here k± = «' =F «' = sinh _1 [sin(^;7/2)] 1 / 2 , k" = 
sin" 1 [sin(^C//2)] 1 / 2 . The uniqueness of this transformation 
requires USt < ir. Incorporating the Trotter decomposition 
and the HS transformation into Eq. (0, the time evolution of 
A is dictated by 



A^O OA 



dsfdsz, ...,dSffI(sf,S2, ...,s%)\.(9) 



The integrand, referred to as as the "Influence Functional" 
(IF), is given by (k = 1, k + p = N) 



1 

2 2 (p+! 



Tr 



K0)a + (4 + p)-^+(4)e iHo(fc - 1)dt e AA e-^( fe - 1 ) 5 ^_(^)...a-(^ +p ) 



(10) 



Here G + (s+) 



,iH St/2 H+{s+) iHoSt/2 



and Q- 



Eq. (O is exact in the St —> limit. Practically, it can be 
evaluated by noting that in standard nonequilibrium situations, 
even at zero temperature, bath correlations die exponentially, 
thus the IF in Eq. (0 can be truncated beyond a memory time 
t c = N s 6t, corresponding to the time beyond which bath cor- 
relations may be controllably ignored 112311 . Here N s is an in- 
teger and the correlation time r c is determined by the nonequi- 
librium situation, roughly r c ~ 1/ Ap. This argument implies 
the following (non-unique) breakup ll23"ll 



I(sf 
xl s (s 



N 



N-N s + V a N-N s +2' 



'2 i •••> °N„ 



( s 2 ' s 3 ' S N B + l) 



It is useful to define the multi-time object 



where each element in the product, besides the first one, is 
given by a ratio between truncated IF, 



-fs(Sfc) Sfc+l, Sk+N s -l) = 



-^( s fe' s fc+l> — ) s k+N s -l) 



k J *fe+H 



'fc+JV 



-2) 



(12) 



'k+2> 



7( S ± 



'2 ' — J °JV 



)-^s( s 2 ' S 3 ' 



■•' s N s + l)-' 



xl s (s 



= 
k 

± „± 



(13) 



and time-evolve it by multiplying it with the subsequent trun- 
cated IF, then summing over the intermediate variables, 



(11) 



fc+2' "fe+3' °k+N s 



E 



± 

fc+i' 



• •) s fc+Ar s -l)^s( s fe+l' s fc+2' — ) S fc+Af s )- 

(14) 

Summation over the internal variables results in the time local 
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expectation value, 

( e *- ^) = ^ 1Z(s k+2 _ N ^,s k+3 _ Ns , ...,s k )(15) 

s ± s ± 

This procedure should be repeated for several (small) values 
of A. Taking the numerical derivative with respect to A, the 
expectation value (A(tk)) is retrieved. 

The main element in this procedure, the truncated IF [Eq. 
(flOt l. is calculated using a fermionic trace formula ll35ll . 

I = Tr[e M ie M \..e M »(p L ®p R ®a(0))] 
= dct{ [I L - f L ] ® [I R - f R ] ® [I s - f s }. 

+ e m ^e m \..e m *[f L ®f R ®fs]}. (16) 

Here, p a , the time-zero density matrix of the a = L, R 
fermion bath and <t(0), the subsystem initial density matrix, 
are assumed to follow an exponential form. Other terms e M , 
with M a quadratic operator, represent further factors in Eq. 
([Tol l. In the determinant, m is a single-particle operator, cor- 
responding to the quadratic operator M = J2i j( m )i-j c l c j> 
c\ (cj) are fermionic creation and annihilation operators, ei- 
ther related to the system or the baths. The matrices I a and 
Is are the identity matrices for the a space and for the sub- 
system, respectively. The functions fi and f R are the bands 
electrons' energy distribution, f a = [e^^-Ma) + with 
the chemical potential p a and inverse temperature /3 a . 

The determinant in Eq. ( TToT ) is evaluated numerically by 
taking into account L s electronic states for each metal. This 
discretization implies a numerical error. However, we have 
found that with L s ~ 100 states we can reach convergence in 
the time interval of interest. Other sources of error, elaborated 
and examined in Refs. I23ll24ll . are the Trotter error, originat- 
ing from the approximate factorization of the total Hamilto- 
nian into the non-commuting Hq (two-body) and Hi (many- 
body) terms, and the memory error, resulting from the trunca- 
tion of the IF. Convergence is verified by demonstrating that 
results are insensitive to the time step and the memory size, 
once the proper memory time is accounted for. 

As we show below, distinct observables may require dif- 
ferent memory time r c for reaching convergence: The dots' 
occupation and the real part of the subsystem off-diagonal 
element, dia 12 , converge for t c ~ l/Ap. In contrast, the 
charge current and 3eri ; 2 require a memory time at least twice 
longer, as these quantities are sensitive to the bias drop at each 
contact, rather than to the overall voltage bias. It is important 
to note that this scaling is approximate, and the actual mem- 
ory time further depends on the subsystem (dots) energetics 
in a complex way: First, the memory time depends on U in a 
nontrivial manner l36Tl . In the absence of U INFPI numerical 
results are exact, irrespective of the memory size used in the 
simulation. This can be seen from Eqs. ( TTOb and dTTb . where a 
cancellation effect takes place leaving free propagation terms 
only, from t = to the current time. At infinitely large U 
one expects again superior convergence behavior, as simulta- 
neous occupancy is forbidden Pol . Second, the position of 



the dot states with respect to the left and right chemical po- 
tentials affect the convergence behavior. We generally found 
that when the dot states are located within the bias window 
a shorter memory time is required for reaching convergence, 
in comparison to the case where the dot energies are out-of- 
resonance with the bias window. This could be rationalized 
by noting that the decorrelation time for electrons within the 
bias window is short relative to the characteristic timescale of 
electrons occupying off -resonance states. 

III. INFPI NUMERICAL RESULTS 

We present here data for the coherence dynamics <7i = 
(d\ (t)d,2 (t)) and the charge current (7 e ) within the interacting 
double-dot AB interferometer. As we show below, we find 
that finite e-e interactions do not destroy the general charac- 
teristics of the coherence behavior, for the cases U/Y < 4 
considered here. We focus on the following set of parame- 
ters: The double-dot subsystem includes two degenerate states 
with e = e m (m = 1,2). The dynamics is studied away from 
the electron-hole symmetric point, E4 = e + U/2 = 0.2. 
The metals' band structure is taken identical at the two ends, 
and we use leads with constant density of states and a sharp 
cutoff at D = ±1. The inter-dot repulsion is taken at the 
range U = — 0.2, whereas the system-bath hybridization 
strength (see definitions in Sec. II) is taken as Y = 0.05. 
As we demonstrate below, our results generally converge for 
U/Y < 4. The bias voltage is applied in a symmetric man- 
ner, pl = —p R , and we take p^ — p R = Ap ~ 0.6. The 
temperature is varied, where /3 = 1/T = 200 corresponds to 
the low-T case, and (3 = 5 reflects a high-T situation. The 
numerical parameters adopted are L s ~ 100 states per bath, 
time step of St ~ 0.8 — 1.6 and a memory time t c ~ 3 — 10. 
This choice of bath states suffices for mimicking a continu- 
ous band structure l23ll24ll . Also, recurrence effects are not 
observed before tY ~ 10. For simulating dynamics beyond 
that time larger reservoirs are constructed, as necessary. The 
time step was selected based on two (contrary) considerations: 
(i) It should be made short enough, for justifying the Trotter 
breakup, StU < 1. (ii) For computational reasons, it should 
be made long enough, to allow coverage of the system mem- 
ory time with few terms, N s < 8, recalling that r c = StN s . 

Before presenting our results we explain the initial con- 
dition adopted. At time t = the double-dot levels are 
both empty, while the (decoupled) reservoirs are separately 
prepared with occupation functions obeying the Fermi-Dirac 
statistics at a given temperature T and bias. 

A. Coherence dynamics at U = 

We begin by presenting results for the noninteracting case, 
U = 0. Figures [2] and [3] display the time evolution of the 
real and imaginary parts of ai^{t), respectively, for relatively 
large bias Ap = 0.6 and at low temperature. We find that 
^ffo - ! 2 decays at a flux dependent rate after the initial rise. The 
imaginary part, displayed in Fig. [3] saturates with a time scale 
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FIG. 2: Time evolution of the states coherence, in the absence of 
electron repulsion effects. Shown is the real part of 01,2 (t), plotted 
for the phases cj> ranging from to 2ir, top to bottom. Ed = 0.2, 
T = 0.05, U = 0, Afi = 0.6, P = 200, L s = 240. 
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FIG. 3: Time evolution of the imaginary part of ai,2(t), in the 
absence of electron repulsion effects. The phase factors <j> range be- 
tween — 7r to 7T, bottom to top. Other parameters are the same as in 

Fig.m 




FIG. 4: Time evolution of a 1<2 for U = (full line), U = 0.1 
(dashed line) and U — 0.2 (dotted line). Main: Real part of ai^it). 
The three top lines were simulated for </> = 0. The bottom lines were 
obtained using <j> = 7r/2. The numerical parameters are dt = 1, 
TVs = 6 and L s = 120. Inset: Imaginary part of 01,3 (i) when 
<j> = 7r. Numerical parameters are <5i = 1.6, N s = 6 and L s = 120. 
Other parameters are E d = 0.2, Y = 0.05, A/x = 0.6 and /3 = 200. 
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T ll2bl_p efining £7i, 2 (t) = |a-i )2 (i)|e^W, it was argued in 
Ref. 12011 that this relative phase localizes to the values <p = 
— ir/2 or 7r/2 in the long time limit when <j> 7^ 2pir, p is an 
integer. This localization behavior is expected only when the 
(degenerate) dot levels are symmetrically placed between the 
chemical potentials, i.e., for e = 0. Away from this symmetric 
point, using e = 0.2, Fig. [2] shows that the real part of ai t 2 
is finite in the asymptotic limit for any phase, besides ir. It 
is interesting to note though that when </> 7^ 2pir, 5ft(Ti.2 still 
approaches a certain-fixed value, irrespective of the magnetic 
flux. 



B. Coherence and current at finite U 

We now investigate the role of e-e repulsion effects on the 
coherence behavior. Fig. [4] displays the real part of 0X2(7;) 
for two phases, = and <fr = tt/2, and its imaginary part 
for 4> = 7r (inset), for three values of U . Data for ^sai^it) 
atU = 0.2 has not yet converged for the r c adopted, see text 
following Fig. \TU\ In comparison to the U = case, we find 
that general trends are maintained, though the long time co- 
herences are larger in the finite U case. Note our convention: 



FIG. 5: Time evolution of the states coherence for U = 0.1, for 
several phase factors. E d = 0.2, F = 0.05, A^ = 0.6, /3 = 200. 
The real part of o\$. was obtained with St = 1 and N s = 6; the 
imaginary part was simulated with St — 1.6 and N 3 = 6. 



the parameter Ed = e+ U/2 is maintained fixed between sim- 
ulations with different values of U. The trajectory simulated 
extends up to Tt = 3, where convergence is satisfactory. Dif- 
ferent memory times were used for simulating the real part of 
the coherence, and its imaginary part. We adopted r c ~ 5 
when simulating SRoi^, whereas r c ~ 10 was used for acquir- 
ing 5(Ji.2. A more detailed discussion of convergence issues 
is given in Sec. III.C. 

Fig. [5]presents 01.2(7;) for several phases </>, at U = 0.1. 
By comparing the data to the zero-[/ case (see Figs. |2]and[3]), 
we conclude that the symmetry of the off-diagonal elements 
is maintained in the presence of U. The general pattern of the 
coherence is displayed in Figs. |6]and|7] plotting the behavior 
of CT12 as a function of the phase factor, at a particular time, 
Tt = 2, for U = 0, 0.1, and 0.2, at different temperatures. It 
should be noted that by this time the real part of the coherence 
has not yet reached its steady-state value. We find that the 
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coherence symmetry around <f> = tt (for 3fJcr 1 _ 2 ) or <f> — (for 
2) is maintained, though the absolute numbers change. 
Interestingly, while the effect of the temperature is significant 
for 3(Ti 2, showing visible reduction in values at high T, the 
real part of 01,2 is only lightly affected by the temperature. 
The downfall of 3o"i,2 with temperature is also reflected in 
the behavior of the charge current, as we show next. 
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form at low temperatures. The electronic states at the right 
lead in the bias window are not fully empty any longer. Simi- 
larly, at the left lead electronic states overlapping with Ed may 
be empty. Overall, this results in the reduction of the current 
at high T. 
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FIG. 6: Effect of finite U on the coherence. Sftcri^ is plotted as 
a function of the phase factor <j> at a particular time, Tt = 2, for 
different U -values and temperatures. Other parameters are the same 
as in Fig. [5] 



FIG. 8: Charge current through an AB interferometer at low tem- 
peratures, p = 200 (left panel) and high temperatures /3 = 5 (right 
panel) for cj> — 0, it/2, and tt, top to bottom with U = 0.2 (full line), 
U — 0.1 (dashed line), (7 = (dotted line). Other parameters are 
the same as in Fig. [5] Numerical parameters are 5t — 1, N s =6 and 
L s =120. 
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FIG. 7: Effect of finite U on the coherence. 3(71,2 is plotted as a 
function of <j> at a particular time, rt = 2, at finite U and for different 
temperatures. Other parameters are the same as in Fig. [5] 

We study the behavior of the charge current at different 
phases, for different e-e repulsion strengths and temperatures. 
Fig. [8] shows, as expected, a destructive interference pattern 
for electron current in the long time limit when cp = tt, irre- 
spective of the value of U. This perfect destructive interfer- 
ence indicates that charge transport is fully coherent in this 
model. The temporal behavior does show however a sensitiv- 
ity to the value of U, manifesting that systems with variable 
U differently respond to the initial condition. 

In the steady-stat limit the current scales like (J e ) oc [1 + 
cos((/>)], for finite U 05|]. This relation does not hold in the 
short time limit. It is interesting to note that irrespective of 
U and the phase factor the current approaches the steady-state 
limit on a relatively short timescale, Tt ~ 2. At high temper- 
atures, Fig. |8]manifests that the system is still fully coherent, 
while temporal oscillations are washed out. The reduction of 
the current at high temperatures can be attributed to the soften- 
ing of the contacts' Fermi functions from the sharp step-like 



C. Convergence analysis 

We exemplify here the convergence behavior of the real 
and imaginary parts of ai.2 at low temperatures, as well as 
the behavior of the current. Fig. [9] demonstrates that SRci^ 
nicely converges for U = 0.2, for t c > 5. The asymptotic 
limit is practically reached, within ~ 1.5% error, already for 
r c ~ 1/A/x. We confirm that the results are insensitive to 
the particular time step selected (inset). We have also verified 
(not shown) that simulations performed with different phase 
factors similarly converge. 

The convergence of 3cri.2 is generally slower, as we show 
in Fig. [10] While Sici^ converges for r c > 1/A/i, we 
find that 9(71,2 requires memory time at least twice longer for 
achieving convergence. For U = 0.1 3<7i.2 is converging. In 
contrast, at stronger interactions, U = 0.2, the large time step 
adopted results in a Trotter error buildup, and the results seem 
to diverge around r c ~ 10 — 12 (inset). 

We also present the behavior of the charge current at differ- 
ent t c values, see Fig. [TT] It generally converges when r c ~ 6, 
irrespective of the phase factor (not shown), for U/T < 4, in 
agreement with earlier studies |24] . 

Overall, we conclude that we can faithfully simulate the 
time evolution of the coherence 01,2 and the current for A/i = 
0.6 and U/T = 2. For larger U, the real part of a\ t 2, the dot 
occupation, and the current can be still converged 1I23I l24ll . 
The simulation of 3<7i,2 requires longer r c and a shorter time 
step at U/T > 2. Roughly, these observations can be ratio- 
nalized noting that the dynamics of SRci ^ is influenced by the 
full potential drop, p^ — pu, similarly to the dots occupation 
(n m ) lHot]- In contrast, the dynamics of 3<7i.2 is sensitive to 
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the bias drop at each contact (2 
lation times. 



, resulting in longer decorre- 
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FIG. 9: Convergence behavior of 5foi,2 for = and U =0.2. Other 
physical parameters are the same as in Fig. [5] Numerical parameters 
are St = 0.8 and N s = 2 (+), N s = 3 (dashed-dotted line), iV s = 4 
(dashed line), N s = 5 (full line) and JV S = 6 (dotted line). The inset 
zooms on the convergence at a particular time, Ft — 2, as a function 
of the memory time r c = N s 5t, using three different values for the 
time steps, St = 0.8 (o), St = 1 (□) St = 1.6 (*). 
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FIG. 10: Convergence behavior of Srcri,2 for (j> = T and (7=0.1. 
Other physical parameters are the same as in Fig. [5] Numerical 
parameters are St = 1.6 and N s = 2 (+), iV s = 3 (dashed-dotted), 
iV s = 4 (dashed line), N s = 5 (full line) and iV s = 6 (dot), N s = 7 
(dotted line). The inset presents Sierra at a particular time, Ft = 2, 
for (7 = 0.1 and U = 0.2, as a function of the memory time r c = 
N s St, using three different time steps, St = 0.8 (o), St — 1 (□) 
St = 1.6 (*). 
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FIG. 11: Convergence behavior of the charge current, <f> = n/2 and 
U =0.2. Other physical parameters are the same as in Fig. [5] Numeri- 
cal parameters are St — 1 and N 3 = 2 (+), N 3 = 3 (dashed-dotted), 
iV s = 4 (dashed line), N s = 5 (full line) and iV s = 6 (dot), N s = 7 
(dotted line). The inset presents the data at a particular time, Ft = 2, 
for f7 = 0.1 (bottom) and U = 0.2 (top), as a function of the mem- 
ory time t c — N s St, using St — 1 (□) and St = 1.6 (*). 



labels the double-dot charge states in order of an empty dot 
(a), single occupied dot, on either the "1" or "2" sites (b and 
c states, respectively) and the state (d), with the two dots oc- 
cupied. Explicitly, \a) <-» |00), \b) <-» 1 10), |c) o 1 01), and 
|d) O 1 11). The creation and annihilation operators of the 
dot are related to this states by d\ <H> |00)(01| + |01)(11| and 
4 O |00)(10| + |10)(11|. Since d\d 2 O |01> <01 1, we iden- 
tify the observable of interest <Ti£=Tr[pd\d2] by a btC - In the 
noninteracting (U = 0) case, the following equations hold in 
in the infinite bias limit |0] 

&a,a = ~^LO- a ,a 

+ 2V R (<T b ,b + <T c ,c + cr b , c e^ /2 + a c . b e" 1 
2T L a a ^ a — 2{F R + T L )a btb + 2T R Gd,d 

2T L a a ^ a — 2(T R + T L )a c _ c + 2V R a d ^ d 
8T*e i<t> f 2 (j b , c + STe'^' 2 ^ 
2T L (a btb + a c , c - e^Vfe 

2r L e^/ 2 a a , a + Sr(a b , 6 + a CjC )e- 
2T R a d4 e~^ 2 - 2(T L + T R )a bt , 



Cc,, 



0~d,d 



0~b,c 



id>/2 



i<j>/2 



(17) 



Here ST = (e^r L - T R ). The hybridization strength, 
independent of the site index m, is defined as T a — 

^Tlik V£ k m^( e ~ e fc)- ^ e equations are valid in the infi- 
IV. MASTER EQUATION ANALYSIS: (7 = AND U = co nhe biafj ^ when ^ ^ | y> r The ^ probabilitV) t0 

occupy any of the four states, is unity, J2j= a bed a j,j = 1- l n 
the steady-state limit we demand that da/dt = 0, the vector a 
includes the matrix elements ak.j of Eq. ( UTt . and obtain the 
stationary solution, valid for <j> ^ 0, 



Rate equations for resonant transport in interacting multi- 
dot structures can be derived based on the microscopic many- 
body Schrodinger equation l37tl . We support INFPI numerical 
simulations with an analytical study of the system's dynam- 
ics, based on such a master equation description. Specifically, 
we adopt the Bloch-type equations derived in Ref. @], for 
the reduced density matrix of the double-dot system in the 
charge state basis Oj.y (t), j = a, b, c, d. Here the index j 



0~b,c(t 



■sin(0/2). 



(18) 



This expression holds in the symmetric setup, = T R , 
for cf> 2irp, p is an integer. One could formally write 
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Ub, c (t) = |crb jC (i)|e^*\ noting that ip equals ±n/2 in the 
steady-state limit. This "phase-localization" behavior was ex- 
plored in Ref. |20t] : The imaginary part of a b , c depends on 
the magnetic phase factor, maximal for <fi = ir with the value 
1/2. The real part is identically zero. The results of Figs. [2] 
and |3] demonstrate the corresponding behavior at finite bias. 
There, the real part is finite, yet small, approaching a fixed 
value. The imaginary part slightly deviates from the predic- 
tion of Eq. ( fT8l due to the finite bias used. One could also get 
hold of the characteristic rates from the dynamical equation, 
by diagonalizing the matrix M in da/dt = Ma. We gather 
five rates, with two phase dependent rates, oc [1 ± cos(^>/2)]. 
For small <j>, the smallest rate is oc [1— cos(0/2)], in agreement 
with lEoll . It can be also proved that in this noninteractingcase 
the steady-state current scales with (7 e ) oc [1 + cos(<^>)] flg]. 

The dynamics of the coherence, attained from the master 
equation ( TPTI i. is displayed in Fig. [12] for 4> = it/2. In the 
long time limit the real part approaches zero; the imaginary 
part reaches isin(7r/4) = 0.354. INFPI results at zero U 
are also included in dotted lines. Deviations of INFPI simula- 
tions from master equation results can be traced down to the 
finite band used within INFPI, in comparison to the infinite- 
flat band assumed in the master equation approach. For finite 
U, we have found that at large bias, A/i = 2D, INFPI data 
basically overlaps with the U = case (not shown) as the 
system basically stands on the symmetric point. 

In the Coulomb blockade regime, for U/T — > oo, a strik- 
ingly different behavior is expected. Starting with the many- 
body Schrodinger equation, one can again derive the sys- 
tem's equations of motion in the large bias limit while ex- 
cluding simultaneous occupancy at both dots, ad,d = 0. 
The following equations of motion are then achieved IH 

(&a,a + &b,b + a cc = 1), 



the eigenvalues of the rate matrix dT~9b , we note that phase de- 
pendent relaxation rates in the Coulomb blockade regime are 
the same as for noninteracting electrons, see also Fig. Q~2] It 
would be interesting to explore this evolution within the INFPI 
approach. However, as we are currently limited to U/T < 4 
values, this would require an algorithmic improvement of the 
INFPI technique. We believe that such an extension could be 
achieved since the Coulomb blockade case should converge 
faster than the intermediate U limit ll36ll . This issue will be 
tackled in our future work. A related switching behavior was 
observed in Ref. where the current, finite in the nonin- 
teracting case for (f> ^ ir, vanishes in the Coulomb blockade 
regime for any phase satisfying </> ^ 2irp. 
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FIG. 12: Master equation analysis: Real and imaginary parts of ab, c 
at 4> = it/ 2, for [7 = and U = oo, obtained by simulating Eq. (17} 
and Eq. dl9t . respectively. Results from INFPI method with (7=0 
are represented by dotted lines, practically overlapping with master 
equation curves. 



&a,a = -4r i cr Q , Q 



2TR(crb,b + 0"c,c + crt> iC e 



-i^/2\ 



&b.b — 2Ti j aa.a 



ZT R (Tb,b — T R (ab, c e l 
2T L a a . a 

ZT R a c , c — T R {ab,cZ 



i<p/2 



cr c ,be 



&c,be 



-i<b/2\ 



a b , c = 2r L e^V Q 

- r fl e~ i0/2 (cr h , b + a c ,c) - 2T R a b ,, 



(19) 



For a spatially symmetric junction, Tl = T R , the steady-state 
solution for the coherence is 

a b , c (i^oo) = -l e ~^ 2 

= -icos(^/2) + |sin(^/2). (20) 

While the imaginary part predicted is identical to the [7 = 
case, see Eq. ( fT8l . the real part is finite and phase depen- 
dent. The Coulomb blockade dynamics is presented in Fig. Q~2] 
(dashed lines). We find that the imaginary part is weakly sen- 
sitive to the onset of U . In contrast, the real part significantly 
deviates from the [7 = case already at Tt ~ 1, By analyzing 



V. SUMMARY 

The intrinsic coherence dynamics in a double quantum dot 
AB interferometer, away from the symmetric point, was sim- 
ulated using an exact numerical technique. At finite inter- 
actions, U/T < 4, at low or high temperatures, we have 
found that the coherence evolves similarly to the [7 = case, 
showing related characteristic timescales and long time val- 
ues. Specifically, we found that for <j> = tt/2 the real part of 
oi,2 approaches a small number (zero at the symmetric point), 
while the imaginary part is larger, ~ 0.35. On the other hand, 
a master equation treatment in the Coulomb blockade regime 
predicts a significantly different behavior: The magnitude of 
SftcTi 2 and 3cri 2 should be the same, V2/4, for the <p = ir/2 
phase factor. 

Future work will be devoted to the study of related models, 
including the spin degree of freedom at each dot. This model 
should demonstrate a decoherence process due to the intrinsic 
spin-flipping dephasing effect Q. Other topics of interest are 
algorithmic improvement of the INFPI technique, to allow for 
the study of the Coulomb blockade regime. One could also 
add a local degree of freedom on one of the interferometer 
arms, e.g., a quantum point contact l38ll or a vibrational mode 
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ll39ll . and observe the time evolution of the interference pattern Award of DS. Fruitful discussions with M. Bandyopadhyay 
in this "which path" experiment. are acknowledged. 
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